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ABSTRACT 

We present UBVRIz' optical images taken with CTI04m+MOSAIC of the 0.32 deg 2 Extended Hubble Deep 
Field South. This is one of four fields comprising the MUSYC survey, which is optimized for the study of galax- 
ies at z = 3, AGN demographics, and Galactic structure. Our methods used for astrometric calibration, weighted 
image combination, and photometric calibration in AB magnitudes are described. We calculate corrected aper- 
ture photometry and its uncertainties and find through tests that these provide a significant improvement upon 
standard techniques. Our photometric catalog of 62968 objects is complete to a total magnitude of R^b = 25, 
with /?-band counts consistent with results from the literature. We select z — 3 Lyman break galaxy (LBG) can- 
didates from their UVR colors and find a sky surface density of 1 .4 arcmin" 2 and an angular correlation function 
w(9) = 2.3 ± 1 .Of? -0 8 , consistent with previous findings that high-redshift Lyman break galaxies reside in mas- 
sive dark matter halos. Our images and catalogs are available at |http : / /www, astro . yale . edu/MUSYC| 

Subject headings: surveys,galaxies:high-redshift,galaxies:photometry 



1. INTRODUCTION 

The study of galaxy formation and evolution requires de- 
tailed information about statistically significant samples of 
dim objects. This, in turn, requires deep imaging and spec- 
troscopy over wide areas of the sky. In pursuit of these 
data, several wide-deep surveys are now underway. Those 
covering several square degrees or mo re either lack spectro- 
scopic follo w-up (e.g. NOAO DWF S. [Tannuzi & DevlH999l 
and ODTS, iMacDonald et all 120041) or are restricted to the 
study of objects at z < 1 (except for quasars) by their imag- 
ing depth e.g. the Slo an Digital Sky Survey (lYork et al.l2000l 
and Abaza iian et alJ 120051)" the VIRMOS-VLT Deep Sur- 
yev (iLe Fevre et al.ll2004t iRadovich et alJl2004l) and DEEP2 
( Davis et al.l I2003|K Other surveys target the high-redshift 
universe with deep HST imaging over fract ions of a square 
degree i.e. the Hubble Deep Fields C Williams et alJ Il996l 
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2000), the Hubble U l tra Deep Field (HUD F), GOODS 
llGiavalisco et all 12004 bickinson etail 12004 . and GEMS 
(Rixetal. 2004). Spectroscopic coverage is feasible over 
these areas but is presently unable to probe deeper than R ~ 
25, making the added imaging depth useful only for morpho- 
logical studies and photometric redshifts. 

The Multiwavelength Survey by Yale-Chile (MUSYC) 
probes the intermediate regime of a square degree of sky to 
the spectroscopic limit of R ~ 25. Section 2 describes the de- 
sign of our survey. §3 reports our imaging observations for the 
extended Hubble Deep Field South. §4 describes our imag- 
ing reduction, and §5 covers our photometric calibration and 
photometry. §6 gives our results for /?-band number counts 
and the sky density and angular clustering of f/V/?-selected 
Lyman break galaxies. §7 concludes. Our analyses assume 
a standard ACDM cosmology with £1,„=0.27 ,Ha = 0.73, and 
flo=70km/s/Mpc. 

2. SURVEY DESIGN 

The Multiwavelength Survey by Yale-Chile (see 
|http : / /www . astro . yale . edu/MUSYC) is designed 
to provide a fair sample of the universe for the study of the 
formation and evolution of galaxies and their central black 
holes. The core of the survey is a deep imaging campaign in 
optical and near-infrared passbands of four carefully selected 
30' x 30' fields. MUSYC is unique for its combination of 
depth and total area, for additional coverage at X-ray, UV, 
mid- and far-infrared wavelengths and for providing the 
U BVRIz'JHK photometry needed for high-quality photomet- 
ric redshifts over a square-degree of sky. The primary goal 
is to study the properties and interrelations of galaxies at a 
single epoch corresponding to redshift ~ 3, using a range 
of selection techniques. We chose to use the UBVRIz' filter 
set in the optical in order to obtain six nearly-independent 
flux measurements with the broadest possible wavelength 
coverage. 
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Lyman break galaxies at z — 3 are selected through their 
dropout in f/-band images combined with blue continuua in 
BVRIz 1 (A > 1216A in the rest-fram e) typical of recent star 
formation JSteidel et alJU996HH999l 1200 3!). Imaging depths 
of U,B,V,R ~ 26 were chosen to detect the LBGs, whose lu- 
minosity function has a characteristic magnitude of m* = 24.5 
in Rab, and to find their Lyman break decrement in the U filter 
via colors {U-V)ab > 1-2. 

Lyman a emitters at z — 3 are selected through additional 
deep narrow-band imaging using a 50A fwhm filter centered 
at 5000A. These objects can be detected in narrow -band imag- 
ing and spectroscopy by their emission lines, allowing us to 
probe to much dimmer continuum magnitudes than possible 
for Lyman break galaxies. 

It has recently become clear that optical selection meth- 
ods do not provide a full census of the galaxy population at 
z ~ 3, as t hey miss objects which are faint in th e rest-frame 
ultraviolet jFranx et alJl2003l iDaddi et alJ l2~004). With this 
in mind, MUS YC has a comprehensive near-infrared imaging 
campaign. The NIR imaging comprises two components: a 
wide survey covering the full square degree and a deep sur- 
vey of the central 10' x 10' of each field. This division be- 
tween deep and wide was chosen because of the 10.5' field- 
of-view of the ISPI near-infrared camera on the CTIO 4m 
telescope. The 5cr point source sensitivities of the wide and 
deep components are K s j,b = 22.0 and K s _ab = 23.3 respec- 
tively. NIR imaging over the full survey area provides a criti- 
cal complement to optical imaging for breaking degeneracies 
in photometric redshifts and modeling star formation histo- 
ries. Deeper JHK S imaging over 10' x 10' sub fields opens up 
an additional win dow into the z ~ 3 universe as the J- K se- 
lection technique (Franx et al. 2003; van Dokkum et al. 2003, 
2004) will be used to find evolved optically-red galaxies at 
2 < z < 4 through their rest-frame Balmer/4000A break. 

Extensive follow-up spectroscopy is being conducted over 
the square degree. A subset of the color-selected Lyman break 
galaxy candidates will turn out to be AGN based on broad- or 
narrow-line emissio n features seen in follow-up spectroscopy 
( Steid elet alJ l2002). Damped Lyman a Absorption systems 
JWolfeetal.l M986) at z > 2.3, which comprise the neutral 
gas r eservoir needed to form most of the stars in the uni- 
verse ( Wolfe et al. 2005), will be searched for in the spectra of 
the brightest color-selected LBG/AGN candidates (typically 
quasars at z ~ 3). 

In addition to the optical and near-infrared, imaging cam- 
paigns at other wavelengths and follow-up spectroscopy are 
integral parts of MUSYC. X-ray selection will be used to 
study AGN demographics ov er the full rang e of accessi- 
ble redshifts, < z < 6, (see iLira et alJ 120041) with Spitzer 
imagin g used to detect optically- an d X-ray-obscured AGN 
( Trei ster et all200HlLacv etaIT20 04). This also allows a cen- 
sus of accreting black holes at z — 3 in the same fields to study 
correlations between black hole accretion and galaxy proper- 
ties at this epoch. 

Future epochs of optical imaging will be used to conduct a 
proper motion survey to find white dwarfs and brown dwarfs 
in order to stud y Galactic structur e and the local Initial Mass 
Function (see Altmann et al. 2004); the additional epochs will 
also enable a variability study of AGN. 

The four survey fields (see Table[TJ were chosen to have ex- 
treme ly low reddening, HI column density ( Burstern" & Heilesl 
119781) . and 100/im dust emission (Schlegel et al. 1998) in or- 
der to facilitate satellite coverage with Spitzer, HST, Chandra, 



and XMM, to take advantage of existing multiwavelength data 
and to enable flexible scheduling of observing time during the 
year. Additionally, each field satisfies all of the following cri- 
teria: minimal bright foreground sources in the optical and ra- 
dio, high Galactic latitude (\b\ > 30) to reduce stellar density, 
and accessibility from observatories located in Chile. The sur- 
vey fields will be a natural choice for future observations with 
ALMA. 

The remainder of this paper describes our optical images 
and catalog of E-HDFS." The techniques used for data re- 
duction and photometry are the same as those used for the 
analysis of the other three fields. Optical imaging from the 
full survey will be reported in E. Gawiser et al. (2005, in 
preparation). The near-infrared data will be discussed in 
R. Quadri et al. (2005, in preparation). The E-HDFS has 
deep public space-based observations at UV, optical, near- 
infrared, and far-infrared wavelengths. The HDFS itself cov- 
ers a small - 2.5' x 2.5' central region with WFPC2 plus STIS 
and NICMOS regions, with deep ground-based JHK cover- 
age of the WFPC2 region available from the FIRES survey 
jLabbeet al J 120031) . Spitzer IRAC and MIPS coverage of 
the central 5' x 15' is being performed in GTO time. The 
extended area aroun d H DFS has previously been imaged by 
iPalunas et all (120001) and lTeolitz et alJ ( 1200 ll) to a depth suf- 
ficient for the study of galaxies at z < 1 . These images were 
made public and combined with deep H -band images by the 
Las Campanas Infrared Survey (LCIRS, Chen et al. 2002) to 
study red galaxies out to z — 1.5. Our survey goes about one 
magnitude deeper in UBVRIz' to probe the z = 3 universe. Our 
Extended Hubble Deep Field South (E-HDFS) field center 
(see Tabled was chosen to keep a bright star (m=6.8) which 
lies just North of the WFPC2 field off of the CTIO+MOS AIC 
detectors. LCIRS covered an H-shape centered on WFPC2 
and thereby provides public H-band coverage of roughly half 
of our E-HDFS field. 

3. OBSERVATIONS 

Optical images of E-HDFS were taken on the nights 
of 2002 October 6,8,10,12 and 2003 May 26,27,28 using 
the 8192x8192 pixel MOSAIC II camera consisting of 8 
2048x4096 CCDs, each with 2 amplifiers, on the Blanco 4m 
telescope at CTIO. Afternoon calibrations were obtained, in- 
cluding zero-second exposures to trace the readout bias pat- 
tern, and domeflats in each filter to be observed except for U 
for which the counts from the dome lamps were insufficient. 
Twilight flats were therefore obtained for the U filter. Dark 
exposures of comparable length to the imaging observations 
were obtained, but due to the negligible dark current they were 
not used for data reduction. A standard dither pattern was 
used to fill in gaps between the 8 CCDs. The pixel scale is 
0.267"/pixel, leading to coverage spanning 37' x 37' of sky 
with each pointing. Figure[0shows the filter response curves 
and their multiplication with the CCD quantum efficiency and 
atmospheric transmission at one airmass. Table[2]gives details 
of the exposure times in each filter for each run along with the 
approximate average seeing measured in the raw images dur- 
ing observations. The seven nights were mostly cloudless, but 
moderate clouds affected some of the V imaging of E-HDFS 
on 2002 October 10,12 and some of the R imaging on 2003 
May 28; our reduction methods described below allow these 
images to be used without biasing the photometry. Photomet- 

2 The data presented here are available at 

|http : / /www .astro . yale ■ edu/MUSYC| 
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FIG. 1.— MOSAIC II UBVRIz' filter set, starting with U at left. Dotted 
lines show filter throughput, and solid lines show total system throughput 
after multiplying by CCD quantum efficiency and atmospheric transmission 
at one airmass. 



ric standard fields from Landolt ( 1992) large enough to cover 
the full MOSAIC II field of view were observed each night, 
and the nights of 2002 October 6 and 2003 May 26 proved 
photometric. 

4. DATA REDUCTION 

UBVRIz' images from CTI04m+MOSAIC II were reduced 
using the MSCRED and MSCDB packages in IRAF v.2.12 3 
following the NOAO Deep Wide Field Survey cookbook 
v.7.02. 4 We used custom software to work around a few diffi- 
culties in these packages, as described below. 

A composite zero image is subtracted from each raw im- 
age to remove the amplifier bias level and pattern. The re- 
sulting image is then "flat-fielded" by dividing by a compos- 
ite domeflat, or a composite twilight flat in the case of U- 
band. A superskyflat for each filter is made by combining 
all of the flat-fielded, unregistered images taken in each filter 
each night with rejection used to remove sources. Our 2' am- 
plitude dither pattern was designed to eliminate the wings of 
bright sources from the superskyflat. Each flat-fielded image 
is then divided by the appropriate superskyflat, which offers 
an estimate of the pixel-by-pixel response to the spectrum of 
the night sky with sufficient counts to achieve 1% precision 
per pixel. Because the superskyflat was produced using flat- 
fielded images, dividing by it serves to remove the original 
domeflat (twilight flat) from the reduction process. The real 
influence of the original flat-fielding is to remove the illumi- 
nation pattern and gross pixel-to-pixel variations before look- 
ing for cosmic rays and bright objects to reject in making the 
superskyflat. Using the superskyflat to correct for the pixel re- 
sponse is considered preferable to using the domeflat because 
the CCD response to the spectrum of the dome lamp (or twi- 
light sky) can have significant systematic differences from its 

3 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 

4 The NDWFS cookbook can be found at 
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response to the spectrum of the night sky. For background- 
limited photometry, this is an important effect. 

We then find an astrometric solution for each image, start- 
ing with fiducial WCS headers provided for each fits exten- 
sion of the raw images and comparing the claimed positions 
with the known positions of stars in the USNO-B catalog us- 
ing the MSCRED routine msccmatch. We found msccmatch 
to be finicky; for some runs the fiducial headers were too inac- 
curate to be corrected with the maximum second-order terms 
used by the MSCRED package. An iterative, non-interactive 
procedure of calling msccmatch multiple times proved suffi- 
cient. The final rms astrometric errors are between 0.2" and 
0.3" in each image, which is consistent with the uncertain- 
ties on individual USNO-B stars despite having fit many more 
stars than free parameters, and the actual solution should be 
better than 0.2" rms. 

To perform aperture photometry later on, our final images 
are transformed to have a common pixel scale and tangent 
plane projection point. This is accomplished by project- 
ing each processed image (after bias-subtraction, flat-fielding, 
superskyflat-fielding, and correction of WCS header informa- 
tion) onto the tangent plane of a common reference image 
using mscimage. The reprojection performed by mscimage 
should be performed with flux conservation set to no the first 
time but then set to yes for any further reprojections. This 
is because the pixel scales in raw MOSAIC II images are 
a function of radius from the field center with several per- 
cent variation from center to corner of the field. The process 
of flat-fielding removes the true illumination pattern caused 
by more photons falling in larger pixels leading to an illu- 
sion of flat background counts despite the variation in pixel 
scales. When mscimage is used for astrometric projection into 
a uniformly sized grid of pixels, flux conservation would re- 
introduce the illumination pattern, thereby causing a variable 
photometric zeropoint across the image. Turning off flux con- 
servation instead produces an average of the values of initial 
pixels neighboring each final pixel location that offsets the er- 
ror introduced by flat-fielding. 

We make three different stacked "final" images for each fil- 
ter. Unweighted versions are made first in order to determine 
the weights for the other two versions. The versions denoted 
by xs were combined using weights optimized for surface 
brightness, which is the method used by the NDWFS cook- 
book. The versions denoted by ps were made using weights 
optimized for point sources; details can be found in the Ap- 
pendix. 

We synthesize a BVR_ps image to use as our photometric 
detection image by adding the B, V, and R ps stacks using 
point-source-optimized weights derived from the signal and 
noise statistics of these stacks. We decided not to add U to 
BVR although it is also very deep because some low-redshift 
objects have different morphologies in U, the seeing is typi- 
cally worse in U, and objects at z — 3 will be nearly invisible 
in U so adding it to BVR would just add noise for these ob- 
jects. Note that for similar reasons objects that drop out in 
B and V are somewhat less likely to be detected in our BVR 
image than in a VR or R image alone. We prefer the BVR ap- 
proach to the creation of a \ 2 image advocated by Szala v etall 
( 1999) because it allows us to measure object morphologi- 
cal parameters directly from the same image used for object 
detection, which is not recommended in the latter approach. 
Moreover, the optimally-weighted BVR combination reduces 
the emphasis on single-filter outliers employed by the \ 2 ap- 
.htpniach and avoids its flaw of treating both negative and posi- 
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FIG. 2.— 33' X 34.5' combined BVR image of E-HDFS used for object 
detection. The mild background feature at the top just left of center is caused 
by the reflection of a magnitude 6.8 star just outside the field. 



tive sky fluctuations as evidence of an object. 

We trim the final BVR _ps image to the maximum size re- 
gion that has nearly uniform signal-to-noise ratio in all three 
input images B, V, and R. For E-HDFS, the trimmed image is 
7395 x 7749 pixels or 33' x 34.5'. Then the IRAF task imalign 
is used to shift and trim the other images to match, giving our 
images the common origin and size required for aperture pho- 
tometry with SExtractor. The trimmed images are normalized 
to an effective exposure time of one second. After the pipeline 
reduction, the background is flat to better than 1 % in all filters. 
The remaining low level , large scale fluctuations were sub- 
tracted using SExtractor (Bertin & Arnouts 1996) with mesh 
size BACK_SIZE=64 and median filter B ACK_FILTERS IZE 
set to 6 mesh units. 

The resulting final BVR image is shown in Figure [2] with 
a magnified image of a small section shown in Fig. [3] New 
FITS headers have been added to these final images indicating 
input information for SExtractor photometry runs and other 
routines, specifically: 

SATUR_LEVEL, the empirically determined saturation level 
in each image, which is usually a factor of a few less than the 
apparent saturation level of the brightest stars, 
SEEING_FWHM, the mode of the seeing for a set of bright 
unsaturated stars, 

MAG_ZEROPOINT, the AB magnitude of an object that 
gives 1 count per second derived from photometric calibra- 
tion as described below, 

FLUX_ZEROPOINT, the flux in /iJy corresponding to 1 
count per second, 

TOT_EXPTIME, the total exposure time that each final image 
represents, and 

GAIN_TOT, the number of photoelectrons represented by 
each count in the final image. 

5. PHOTOMETRY 
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FIG. 3. — 2' X 2' zoom of BVR detection image. The seeing in the image 
is 0.95" fwhm. 



5.1. Photometric Calibration 

Our calibration scheme is to take ^5 minute exposures of 
our fields (referred to as "calibration images") in all optical 
filters on photometric nights. This was achieved for E-HDFS 
on Oct. 6, 2002. Photometric calibration is performed using 
a version of the Landolt catalog with a ll magnitudes and col- 
ors converted into the AB95 system of Fukuaita et alJ (119961) 
(referred to throughout this paper as AB magnitudes). We in- 
cluded z' magnitudes for Landolt standard stars in our catalog 
of standards using the formula 

z / ab = /ab-0.59(/?-/) A b+0.16(V-/0ab-0.04 , (1) 

which was generated from the formulae of iFukugita et al.l 
( 1996) and has been corrected for a 0.02 magnitude bias found 
when comparing this prediction wit h measurements o f z' mag- 
nitudes of bright Landolt stars by Smit h et al.N2002l) . The re- 
sulting 0.04 magnitude rms error was used to predict errors in 
these estimates. This procedure offers a significant reduction 
in observing time and a tremendous increase in the number of 
calibration stars versus the traditional method of using spec- 
trophotometric standard stars for filters outside the Johnson- 
Cousins system. Only a few s tars per Landolt fie ld have been 
turned into SDSS calibrators ( Smi th et ai]l2002l) and they are 
typically too bright for the 4m telescope to use without defo- 
cusing. 

Table [3] lists the airmass and seeing for each of our cali- 
bration images along with our photometric solution for the 
magntiude in terms of measured counts per second, 

m = — 2.5log l0 (counts/s) + Z— cX — Colorterm , (2) 

where Z is the magnitude zeropoint, c is the airmass coeffi- 
cient, X is airmass, and the color term is listed in Table [3] 
Our images of Landolt standard fields occupied a range of 
airmasses bracketing the airmasses of these calibration im- 
ages but insufficient to break the degeneracy between zero- 
point and airmass coefficient. We therefore used fixed air- 
mass coefficients, with the resulting uncertainty reflected in 
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the uncertainty in the zeropoint fit, which is <1.5% in each 
filter. Since AB95 was carefully calibrated for the standard 
Johnson-Cousins UBVRI filter set, the conversion of the Lan- 
dolt catalog is quite accurate. However, the filter system used 
at CTIO is not a precise match to Johnson-Cousins, mak- 
ing the use of a color term helpful in photometric calibra- 
tion. The relatively small color coefficients shown in Table 
[^minimize the scatter in our photometric solution by provid- 
ing a better estimate of the true AB magnitude of each Lan- 
dolt star in our observed filter system. We set the AB color 
coefficient to zero when determining photometry for our fi- 
nal object catalog; this places the fluxes on our observed filter 
system rather than standard Johnson-Cousins filters. A color 
correction can be calculated later for applications that require 
fluxes extrapolated to Johnson-Cousins colors. Because pho- 
tometric redshift codes multiply spectral templates by the at- 
mospheric transmission, filter transmission, and CCD quan- 
tum efficiency used for the observations, they already account 
for the origin of the color term and should be given fluxes in 
the observed filter system. The presence of non-power-law 
features such as Lyman and Balmer breaks in these template 
spectra (and real objects) makes it desirable to avoid color 
transformations calibrated with stars. 

We then used the calibration images to create standard stars 
in each filter for E-HDFS and used these new standard stars 
to calculate the zeropoint in our final images. This avoids 
attempting to explicitly correct the photometry for the air- 
mass and cirrus extinction in each image which contributed 
to our final images, allowing us to use images taken in non- 
photometric conditions as part of the final stack. For E-HDFS 
we used a set of 10 stars, each from a different original am- 
plifier, chosen to be bright enough to provide good signal in 
U and yet to avoid saturation in R, I, and zf. Care was taken 
to choose some stars near the field center and some near the 
sides and corners of the field to check for systematic effects. 
The photometry of the new standards was determined using 
the IRAF routine phot. We used the same 14" diameter aper- 
ture for photometry of these standard stars in order to mini- 
mize aperture losses. We repeated this process by using phot 
on our final science images for the same stars. We checked 
for systematic patterns in the offsets but found none, and the 
rms scatter was 0.03 magnitudes, implying a smaller error in 
the mean offset. 

The total exposure time, seeing, photometric zeropoints, 
and 5 a point source detection depths (AB) of our final im- 
ages are given in Table|4] Flux zeropoint s are also determined 
using the definition ( Fukugit aet all 19961) 

fApiy) = 3.631 x 10 9 • 10-° 4mAB . (3) 

Given the measured uncertainties in the initial zeropoints and 
our measured scatter in the solution for the zeropoints of the 
final images, we estimate a total systematic uncertainty of 3% 
in the zeropoint of each filter. 

5.2. Correlated Noise on Small and Large Scales 

Derivations of optimal photometric techniques and their un- 
certainties typically assume Poisson noise which is uncorre- 
lated between image pixels. However, the astrometric repro- 
jection used to place all images on a common grid correlates 
neighboring pixels, leaving the large-scale noise properties 
unchanged but making a pixel-by-pixel rms a poor measure 
of the typical noise due to sky fluctuations in a larger aper- 
ture. As this is the method used by SExtractor to estimate 
photometric errors, we undertook a detailed investigation to 



test its accuracy, which turns out to be poor, and to find an 
improved method. Our empirical method for estimating pho- 
tometric uncertainties as a function of aperture size is also 
sensitive to large-scale noise correlations from unsubtracted 
nearby objects and from any CCD defects that survive flat- 
fielding. Wings from bright objects appear to affect the noise 
statistics; when we set all pixels belonging to detected objects 
to zero and displayed the resultant image, wings are clearly 
visible around the brightest objects. 

We estimate the error due to background fluctuations and 
noise correlations in each filter via a custom IDL code which 
places ~ 2000 random apertures of a given size on the sky- 
subtracted image which do not overlap with any of the pixels 
in the segmentation map of isophotal object regions produced 
by SExtractor. We use circular apertures of area n P j X centered 
at integer pixels and describe them by an effective size N = 
Jn~^~ x . A Gaussian is then fit to the histogram of aperture 
fluxes to yield the rms background fluctuation as shown in 
Fig. |4] The histograms appear well described by the best-fit 
Gaussians. 
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FIG. 4. — Histogram of aperture counts for two different aperture sizes 
showing Gaussian fits. The larger aperture has more variation in counts and 
hence is fit by a Gaussian with larger fwhm and lower mode. 

Figure |5] shows the rms of aperture fluxes vs. N for raw, 
sky-subtracted raw, and final V images. The solid curve gives 
our fit to the noise properties of the final V image using the 
function 

a N = a l aN f3 , (4) 

where a = 0.77, (3=1 .30 for the measured rms pi xel noise of 
(7i=0.0 14. An alternative formula suggested by Lab be et alJ 
J200l . 

(T N = aiaN+a l bN 2 , (5) 

yields an equally good fit, with a = 1 .007, b = 0.05 1 . This for- 
mula shows an explicit sum between contributions from Pois- 
son noise, which is independent from one pixel to another 
and hence has rms proportional to N, and correlated noise 
from fluctuations in background level on scales larger than the 
aperture which yields an rms proportional to A^ 2 . These sim- 
ple, extreme cases are shown as dotted lines starting from the 
value of (7 1 in Figure[5]and are seen to bracket the true behav- 
ior except at very small aperture sizes. Apertures of just a few 
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FIG. 5. — V-band rms fluctuation in counts/s as a function of linear aper- 
ture size N in pixels for a raw image (squares), sky-subtracted raw image 
(diamonds) and the final image (circles). Filled symbols show <ri for each 
image which is calculated pixel-by-pixel rather than using a circular aper- 
ture. For the final image, the solid line gives the fit of Eq. |4]and the dotted 
lines show the limiting cases of rms proportional to N 2 and N. 



pixels in area are affected by the small-scale noise correlations 
introduced by re-pixelization performed during re-projection; 
this effect should actually depress u\ since each final pixel is 
the average over several nominally independent input pixels. 
We prefer the simplicity of Eq. |4]which reflects the reality that 
noise exists on a range of scales leading to an effective power- 
law behavior intermediate between these two extremes. Table 
|5]lists our measured rms-per-pixel o\ and the fit coefficients 
a and (3 for all of our filters; the uncertainties in a and (3 are 
highly correlated but both have been determined to 5% preci- 
sion, whereas the uncertainties in a\ measurements are at the 
1% level. 

A similar analysis of our raw images shown in Figure|5]re- 
veals that the noise in the raw images (squares) is fully corre- 
lated due to the domination of background fluctuations. Once 
a sky-subtraction has been performed (diamonds), the noise 
properties resemble those of the final image except for be- 
ing noisier by a factor roughly equal to the square root of 
the number of exposures. While some of this improvement 
comes from flat-fielding, it appears that the correlated noise is 
due to background fluctuations from objects whose subtrac- 
tion improves with the square root of observing time. This 
is inconsistent with confusion noise from undetected sources, 
which is not expected to be a significant factor even in op- 
tical images this deep. One possibility for the origin of this 
correlated noise is that the global sky subtraction performed 
by SExtractor is insufficient. However, when we mimic the 
BACKPHOTO_TYPE LOCAL mode of SExtractor by esti- 
mating the local sky in a square annulus around each aper- 
ture and subtracting this off we find a slight increase in the 
noise correlations. Hence it appears that the background con- 
tributions from sky gradients and nearby objects have been 
subtracted about as well as one could expect, and that the re- 
maining contributions increase the noise significantly above 
the extrapolation from u\ one would assume for the case of 
uncorrected noise. 



5.3. Optimal Apertures for Photometry 

We wish to determine accurate photometry for a wide 
range of object sizes and brightnesses, with particular atten- 
tion to high -redshift galaxies with half-light radii 0.2-0.3" 
( Steid el et al.ll9 96a)) which are dim and typically unresolved 
in our ~ 1 " seein g and can th erefore be approximated as point 
sources. Indeed, Smail et al. ( 1995) showed that most galax- 
ies at m ~ 25 have half-light radii < 0.3" regardless of red- 
shift. In traditional aperture photometry, the flux of an ob- 
ject is measured by giving constant weight to one region of 
the image and zero weight to the rest, with fractional weights 
used for pixels which fall only partially in the chosen aper- 
ture. The simplest case is a circular aperture centered at the 
object centroid. For the case of uncorrected noise dominated 
by the sky background and a Gaussian PSF, the circular aper- 
ture with optimal signal-to-noise is easy to derive and has a 
diameter equal to 1.35 times the seeing fwhm. However, the 
actual PSF in each filter has a near-Gaussian core with broad 
wings i.e. more of the flux is found outside the fwhm than 
for a Gaussian. To investigate the effect of the non-Gaussian 
PSF on photometry we measured the enclosed flux as a func- 
tion of radius for a set of objects selected to be bright but 
unsaturated (19 < m < 20) and unresolved according to SEx- 
tractor (stellarity > 0.8, as discussed below). Fig. [6] shows 
the results for the U, V, and / bands plus a Gaussian of 0.96" 
fwhm that represents an excellent fit to the core of the V-band 
PSF but only encloses 80% of the total flux. A 14" diameter 
aperture was used to define the total object flux, with rapid 
convergence seen in the bluer filters but only gradually in R,I, 
and z 1 ■ While this could be due in part to contributions from 
neighboring objects, we used the median statistic to define 
these profiles to avoid sensitivity to the small fraction of ob- 
jects with bright neighbors. The contribution of nearby ob- 
jects should lead to a divergence with increasing aperture di- 
ameter which is not seen, so it appears that the PSF wings do 
contain a higher fraction of the total flux in the redder bands, 
with perhaps a few percent missed by our 14" diameter aper- 
ture. 

Figure [7] plots the enclosed flux, the rms background fluc- 
tuation, and the resulting signal-to-noise for point sources as 
a function of aperture diameter. There are competing effects 
which serve to make the optimal aperture slightly smaller than 
the idealized result of 1.35 x fwhm. The non-Gaussian PSF 
places more flux outside the "optimal" diameter, so a larger 
aperture would find more signal. However, the large-scale 
noise correlations mean that noise increases with aperture size 
faster than in the ideal Poisson case. Table[6]lists various see- 
ing measurements along with the optimal circular aperture di- 
ameters and fractional flux enclosed within these apertures for 
each filter which is used to correct the observed fluxes to esti- 
mate the total object flux in each filter. A systematic increase 
in the seeing occurs when the estimate is made using SExtrac- 
tor's fwhm parameter, which assumes a Gaussian PSF, instead 
of IRAF's imexam routine, which allows a Moffat fit. The 
non-Gaussian PSF shape is reflected in the further increase 
in the diameter containing half of the object light implied by 
the half-light radii ru2 output by SExtractor; a Gaussian con- 
tains half of the flux within one fwhm. The optimal circular 
apertures for point sources turn out to be quite close to these 
values for the half-light radii. The image quality is measured 
to be nearly constant across the 33' x 34.5' field, with system- 
atic variations in fwhm and half-light radius seen at the 10% 
level in the / band but at no more than the 5% level in the 
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FIG. 6. — Left panel: Median fractional flux enclosed for bright, unsaturated point sources as a function of aperture diameter. Right panel: PSF implied by this 
enclosed flux profile. In both panels, the upper solid curve shows results for the V band, with the lower solid curve representing a Gaussian PSF which matches 
the median fwhm of 0.96" measured by SExtractor on the V image but only contains 80% of the total object flux. This shows that the PSF core is Gaussian but 
the wings are broader. The long-dashed curve shows results for the U band, which has worse seeing but equally rapid convergence to the total object flux in the 
left panel. The short-dashed curve shows results for the / band, which has nearly identical seeing to V but the slower convergence to the total object flux seen in 
the Rh' filters. 
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FIG. 7. — The long-dashed curve shows the V-band noise as a function of 
aperture diameter, normalized to the value at 14". The short-dashed curve 
shows the V-band enclosed flux from Fig. |S] The solid curve shows the 
resulting signal-to-noise ratio, which has a maximum at 0.5" where less than 
half of the flux is enclosed. The behavior for other filters is similar. 



other filters. 

For extended sources, our apertures are no longer optimal, 
and the corrections for the fraction of enclosed flux are in- 
correct. These sources are barely affected by seeing smaller 
than their intrinsic size, leading to a flux better described by 
the aperture area than by the fractional enclosed point-source 
flux. The traditional solution is to use larger apertures (> 2" 
diameter) to measure colors of extended objects as well as 
point sources. Our analysis shows that the signal-to-noise of 
point source photometry is reduced by up to 30% by using 



these larger apertures, with a bias in the U — V color of point 
sources of 20% (0.2 magnitudes) introduced by the different 
seeing in these images. This could be solved by convolving all 
of the image s to the relatively p oor seeing (fwhm=1.3") of the 
U band (see Labbe et al. 2003). However, the non-Gaussian 
PSF shape requires a non-Gaussian convolution kernel; sim- 
ply matching the fwhm by Gaussian convolution is not suffi- 
cient to fully match the PSF. Moreover, smoothing the other 
images to match the PSF shape of the U band would be costly 
in terms of decreasing the signal-to-noise for point source 
photometry, especially given the observed large-scale correla- 
tions in background noise. In choosing the "optimal" aperture 
sizes for each filter in Table [6] we have chosen amongst the 
range of apertures with signal-to-noise>95% of the optimal 
value in a manner that reduces the range of aperture sizes in 
order to reduce the errors introduced for extended sources. 
This has the effect of choosing slightly larger apertures in 
the redder filters than would formally maximize the signal- 
to-noise. This provides additional benefits, as the approxima- 
tions that high-redshift galaxies are point sources, image co- 
registration is perfect, and object positions have no error all 
become less valid as the apertures decrease in size. Our sim- 
ulations indicate that these combined effects will not affect 
object colors beyond the 10% level for the chosen aperture 
sizes. 

Our primary solution to the challenge of obtaining decent 
colors for both unresolved and extended objects dimmer than 
the night sky is to create corrected aperture fluxes (hereafter 
referred to as APCORR). We tried using the A_IMAGE and 
B_IMAGE parameters to model the flux distribution as an el- 
lipsoid but encountered too many systematic problems with 
these parameters to prefer this approach. We therefore use the 
half-light radius (FLUXJRADIUS) of each object determined 
by SExtractor in the BVR detection image and assume a two- 
dimensional gaussian light profile with twice this value as 
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fwhm. The half-light radius measured in the BVR detection 
image represents the object's intrinsic size convolved with the 
BVR seeing, and this allows us to predict the observed size in 
each filter given the known seeing measured from the median 
half-light radii of bright point sources listed in Table [6] For 
object i in filter j we expect 



2 1 

ij 2 In 2 



2 

ps,j 



'ps,BVR 



>] 



(6) 



where the last term in parentheses modifies the quadrature 
sum of intrinsic object size and BVR seeing to the needed 
sum of intrinsic object size and seeing in filter j and the pref- 
actor converts from half-light radius (which equals half-width 
half maximum for a gaussian) to the rms of the light profile. 
The fractional flux of object i falling within the optimal cir- 
cular aperture radius r ap j for filter j listed in Table|6]is then 
assumed to be 



fracij = min frac ps j , 1 — exp j 



ap,J 



2<rf, 



(7) 



where the first value in parentheses refers to the fractional 
flux for point sources listed in Table [6] and the second gives 
the fractional flux of the two-dimensional gaussian which ap- 
proximates the object. Taking the minimum of these two val- 
ues prevents noisy data from causing an object to be inferred 
to be smaller than a point source. Each object's flux measured 
within the optimal aperture for a given filter is then corrected 
by dividing by fracjj to infer the total flux in that filter. 

Another solution for attempting to handle both extended 
and unresolved sources is SExtractor's AUTO aperture, which 
uses a Kron ellipse adapted to each object's light profile to 
find ~ 94% of the flux for extended sources and ~ 97% for 
point s ources using the parame ters we have adopted, accord- 
ing to lBertin & Arnoutsfj 19961) . In the U, B, and V filters, we 
find that AUTO indeed measures 96% of the point source flux 
found in a 14" diameter aperture, with the value decreasing to 
93% in R, 92% in /, and 89% in z' in keeping with the broader 
PSF discussed above. 

5.4. Estimation of Photometric Uncertainties 

Because SExtractor assumes Poisson sky noise, it produces 
underestimated photometric uncertainties for objects dimmer 
than the background sky. 6 The standard formula used by SEx- 
tractor to predict photometric errors in flux as a function of 
aperture size is 



(7 



phot ,SE 



lpix+ GAIN 



(8) 



where the first term represents sky fluctuations assuming un- 
correlated noise between pixels, F is the object flux in counts 
within the aperture, and GAIN is the total effective gain used 
to turn counts in the normalized image into original photons 
detected by the CCD, with the fraction F / GAIN giving the 
Poisson variance. SExtractor measures a\ using the median 
fluctuation in its internally computed background map, and 
this measurement agrees with ours to within a few percent. 

5 We found FLUX_RADIUS to be considerably more robust than FWHM, 
making it the best measure of the object light distributions. 

6 We found it necessary to run SExtractor with WEIGHT_TYPE BACK- 
GROUND,BACKGROUND in order to get it to use the background in the 
BVR detection image as a source of uncertainty in object detection and the 
background in the measurement image as a source of uncertainty in the pho- 
tometry. 



We modified this to make use of our fits for (correlated) back- 
ground noise which yields 



° 2 phot = ° 2 i<x 2 np ix + -^f N - 



(9) 



where the first term now represents empirical sky fluctuations 
which include the subdominant contribution from electronic 
readout noise. Both noise terms will be larger in regions of the 
image where reduced exposure time leads to a lower effective 
gain and hence higher background rms. SExtractor accounts 
for this by using a variance map of; which is either input as 
a WEIGHT_IMAGE or generated internally from a map of 
background fluctuations in the image. The errors produced 
by SExtractor as a function of position are then 



<-^'{ n P" + a 2. GAlN ) ■ 



(10) 



The ratio of the formulas in equations [8] and [9] allows us to 
make a simple correction to the flux uncertainties output by 
SExtractor to yield 



a 2 nL+- F 



P« a, -GAIN 



a 



phot ~ a phot.SE 



n„ix + 



(ID 



P' x a 2 , -GAIN 



For negative fluxes, F is set to zero in all of these error terms, 
since the uncertainty is entirely dominated by the background 
fluctuations. For AUTO apertures, the area is given by that of 
the Kron ellipse used i.e. n pix = 7rr^ on A_IMAGE B_IMAGE 
(the so-called KRON_RADIUS is really a scaling of the 
A_IMAGE and B_IMAGE parameters based on the second 
moment of the object light distribution). 

The resulting rms has units of counts per second and is con- 
verted to flux using the zeropoints in Table|4] The uncertainty 
determined for optimized aperture fluxes is then multiplied 
by the same correction factor used to estimate the total object 
flux i.e. the reciprocal of frac,j defined in Eq. Finally, 
the uncertainty in the corrected aperture flux is increased by 
a factor (frac ps j/fracij) 2 which serves to amplify the uncer- 
tainties for extended sources to account for the uncertainty in 
their correction factors. This reproduces the errors found in 
corrected aperture flux es for extended objects in our simula- 
tions described in 35.51 The photometric errors derived in this 
manner do not include the 3% calibration uncertainty which 
is common to all sources in a given filter. 

5.5. Photometry Tests on Simulated Sources 

We used the IRAF package artdata to simulate stars and 
galaxies with known magnitudes and positions and to add 
them to our observed images to get realistic crowding effects 
and background noise. 7 We ran SExtractor on these "simu- 
lated" images in dual-image mode to simulate our full pho- 
tometric pipeline. We found object centroiding errors to be 
< 0.3 pixels i.e. < 0.1" at R < 25. in the centroiding of indi- 
vidual objects We found AUTO photometry to be nearly unbi- 
ased for both single filter fluxes and colors. However, Figure 
[8] shows that AUTO has larger errors for point sources than 
the corrected aperture fluxes due to the larger AUTO aper- 
tures including significantly more sky noise, so we recom- 
mend AUTO fluxes only for extended sources and give fluxes 

7 The number of simulated sources was roughly one-tenth that of real 
sources so the crowding characteristics of the images can be considered un- 
changed. 
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FIG. 8. — Errors, i.e. (observed - true), in estimated V fluxes versus errors in estimated U fluxes for simulated point sources. Left panel shows results for 
the AUTO aperture as chosen by SExtractor, and the right panel shows results for the corrected aperture fluxes described in the text. The APCORR fluxes are 
unbiased and have a significantly smaller scatter than AUTO. 
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FIG. 9. — Errors in estimated V fluxes versus errors in estimated U fluxes for simulated galaxies. The left panel shows results for the AUTO aperture as chosen 
by SExtractor. The right panel shows results for the corrected aperture fluxes described in the text. APCORR appears nearly unbiased, and the size of the errors 
for galaxies are comparable to those of AUTO but exhibit different systematics. 



for both types of apertures in our catalog. The errors seen in 
the U and V filters appear to be entirely uncorrelated. Figure 
[9] shows that AUTO performs similarly to corrected aperture 
fluxes for extended objects, although both show significant 
covariance between U and V due to misestimation of the true 
object light profile. 

Figure [10| shows the flux errors divided by their uncertain- 
ties for simulated point sources. APCORR again appears un- 
biased but this plot reveals that the error estimates are roughly 
equally accurate for both cases. The median squared value is 
close to one in both cases, so the reported flux uncertainties 
in the catalog appear trustworthy for point sources. FigurefTTI 
shows the flux errors divided by their nominal uncertainties 
for simulated galaxies. This plot illustrates that the estimated 



uncertainties are typically too small for galaxies in both cases. 

Figurefl2lshows that U — V colors of simulated objects ap- 
pear slightly biased for AUTO but unbiased for APCORR, 
and the errors are a bit smaller for APCORR for point sources 
and of similar size for both types of fluxes for galaxies. The 
color errors are highly correlated, especially for galaxies, im- 
plying that these errors are primarily caused by errors in the 
isophotal object detection by SExtractor on which both fluxes 
depend. These include problems caused by blending with 
nearby neighbors, although the median object is unaffected by 
neighbors in this uncrowded field. Th e AU TO and APCORR 
colors of real objects are discussed in 36.21 

5.6. Aperture Photometry 
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FIG. 10. — Errors, i.e. (observed - true), in estimated V fluxes divided by reported lcr flux uncertainties in V vs. errors in estimated U fluxes divided by reported 
lcr flux uncertainties in U for simulated point sources. Left panel shows results for the AUTO aperture as chosen by SExtractor, and the right panel shows results 
for the APCORR fluxes described in the text. 
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FIG. 1 1 . — Errors in estimated V fluxes divided by reported lcr flux uncertainties in V vs. errors in estimated U fluxes divided by reported lc flux uncertainties 
in U for simulated galaxies. Left panel shows results for the AUTO aperture as chosen by SExtractor, and the right panel shows results for the APCORR fluxes 
described in the text. 



We performed aperture photometry using SExtractor in 
dual-image mode with BVR as the detection image and each 
final image as the measurement image, leading to a common 
catalog for all filters. Detection was performed after filtering 
with a 7 x 7 pixel gaussian convolution kernel with 4 pixel 
( 1 .07") fwhm which is well-matched to the seeing in our BVR 
image (see Appendix). We used a threshold of 1.2 for both 
detection and analysis; a single pixel at the 5 a level in the 
original BVR image would barely make the 1 2a threshold af- 
ter filtering. Using the identical parameters, we detected 318 
objects in a negative of the BVR image, so given the symmet- 
rical nature of sky fluctuations and readnoise we estimate an 
equivalent number of spurious objects in our final catalog i.e. 
0.5%. For each filter, we used the "optimal" apertures and 



corrections described above as well as the AUTO aperture de- 
termined for each object by SExtractor centered on the BVR 
barycenter. The photometr y was corrected for Galactic dust 
extinction of E(B-V)=0.03 (ISchlegel et al.ll998l) . 

The 5 a point source detection depths in Table [6] were de- 
termined by multiplying the a of the measured sky noise for 
the optimal point source aperture sizes by 5 and then correct- 
ing for the flux missed by these apertures to turn sky fluc- 
tuations into total point source magnitudes. Many estimates 
of point source detection limits in the literature are based on 
extrapolating the rms pixel noise, u\ , to the chosen aperture 
size assuming uncorrelated sky noise. For correlated noise 
like that in our images, which we expect is typical, this will 
overestimate the true depth in a 2" diameter aperture by 0.4 
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FIG. 12. — Errors in U-V colors (AB magnitudes) for AUTO aperture as chosen by SExtractor vs. U-V color errors (AB magnitudes) for the APCORR fluxes 
described in the text. Left panel shows results for simulated stars, and right panel shows results for simulated galaxies. 



magnitudes and in a 3" diameter aperture by 0.6 magnitudes. 

5.7. The Photometric Catalog 

Our optical photometric catalog of 62968 objects in 
E-HDFS is available in full in the electronic version of the 
journal and from our website, with the first several lines 
shown in Table [7] All positions and shape parameters are 
measured from the final BVR image. The columns in the 
catalog offer the following information: 
Column 1: Object number starting with first object as 
Column 2: Object name (MUSYC-) 

Column 3: Right ascension in decimal hours, double preci- 
sion 

Column 4: Declination in decimal degrees, double precision 
Column 5: x barycenter 
Column 6: y barycenter 
Column 7: stellarity classification 
Column 8: half-light radius 

Column 9: rms of flux distribution along major axis 
Column 10: rms of flux distribution along minor axis 
Column 11: object position angle counterclockwise from 
x-axis 

Column 12: flags output by SExtractor 
Column 13,15,17,19,21,23: AUTO flux in U,B,V,R,I,z' 
Column 14,16,18,20,22,24: AUTO flux uncertainty in 
U,B,V,R,I,z' from iTTl 

Column 25,27,29,31,33,35: APCORR flux in U,B,V,RJ,z' 
after division by correction factor based on half-light radius 
described in Eq. [7J 

Column 26,28,30,32,34,36: APCORR flux uncertainty in 
U,B,V,R,I,z' from (II It after division by correction factors 
based on half-light radius described in Eq. [71 and multiplica- 
tion by additional factor for extended sources described above 

Photometric measures in the catalog are given in units of 
flux density in /iJy (1 /iJy = 10~ 32 W m" 2 Hz" 1 ). The conver- 
sion to AB magnitudes is simple using the formula 

AB = 23.90-2.51og 10 / v (12) 



jFukueita et aljfl996l) . This avoids the loss of information 
that comes from SExtractor representing objects with nega- 
tive aperture fluxes as having m = 99 and the confusion that 
results from flux errors consistent with a non-detection being 
turned into extremely large magnitudes. Photometric errors 
are nearly symmetrical in flux but not in magnitudes, even 
for objects with low signal-to-noise, making it acceptable to 
represent the uncertainties with a single number. This makes 
flux density the natural unit to be used for calculating pho- 
tometric redshifts if one wishes to use a \ 2 function to com- 
pute the likelihood as is done by the publicly available codes. 
For color selection, one can either turn color criteria into de- 
sired flux ratios between filters or turn the measured fluxes 
into magnitudes and subtract to obtain colors. 

SExtractor offers a neural network classification of objects 
into "stars" and galaxies; stellarity > 0.8 implies an unre- 
solved profile consistent with the PSF, stellarity < 0.3 im- 
plies an extended profile, and 0.3 < stellarity < 0.8 repre- 
sents uncertain objects which are typically too dim for this 
classification to work successfully. Figure IT3lshows the stel- 
larity results for our catalog. The classification starts to break 
down at m > 23, with brighter objects confidently split be- 
tween stellar and extended profiles but many dimmer objects 
receiving an inconclusive value. Stellarity becomes nearly 
useless at m > 25 where the majority of values fall in the un- 
certain regime, which is appro priate given that r l / 2 < 0.3" for 
most galaxies at m > 25 ( Smail et al. 1995). However, the vast 
majority of objects at m > 23 are galaxies due to the rapidly 
rising galaxy luminosity function and the plateau in the stel- 
lar equivalent. This means that for many scientific goals all 
objects with stellarity < 0.8 may be considered galaxies. 

Figure[2]shows the locations of all objects in our catalog. 
There is an overabundance of objects along the image bor- 
der which is not apparent to the eye but can be eliminated by 
removing the 574 objects with flags > 8, which identifies ob- 
jects truncated by the image border. The one visual blemish 
is a line of spurious objects detected along a heavily saturated 
column extending from the brightest star in the field. This line 
is removed by eliminating all 1226 objects with flags > 4, 
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FIG. 13. — SExtractor stellarity classification for objects in our catalog. The left panel shows stellarity versus total R-band AB magnitude and reflects confusion 
for saturated objects at R < 19. The right panel shows a histogram of stellarity, with vertical solid lines separating galaxies (stellarity <0.3), uncertain objects 
(0.3< stellarity <0.8), and stars (stellarity >0.8). The classification is clean at R < 23 (light histogram), shows an increasing fraction of uncertain objects at 
23 < if < 25 (medium histogram), and breaks down at R > 25 (dark histogram). 
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minimized with negligible loss of galaxy counts by requiring 
stellarity < 0.8. The total /?-band magnitude counts (from 
AUTO) are shown in Figure \^\ and appear 90% complete to 
R = 24.5 and 50% complete to R = 25.5. Part of the incom- 
pleteness comes from AUTO begin ning to under estimate the 
flux of very dim objects, as noted by Labbe et al. (2003J). The 
fit for the number of galaxies per magnitude per square degree 
is log(A0 = -3.52+0.34fi. Our slope of 0.34±0.01 agrees well 
with p revious measureme nts of 0.321 ±0.00 1 bv lSmail et alJ 
( 1995), 0.361±0.004 by Capa ket al.1 J200l. 0.31±0.01 and 
0.34±0.01 in two different fields bv LSteidel & Ham ilton 
(11993b . and 0.39 bv lTvsonl 17988). Our differential counts 
atR = 25.5 of 1 .3 x 10 s agree well with the values of 1.2xl0 5 
reported bv lSteidel & Hamilton! lll993l) . 

6.2. UVR Photometric Selection of Lyman Break Galaxies 

Various filter sets have yielded success at finding 
LBGs, i ncluding U n GK iSteidel et al] Il996bl) SPSS 
u'e'r'i 'z' JBentz et all l2004l) . u'BVRI llCooke etaD 12 005 ), 
GKi JSteidel et alJ 11999ft . and BRI (iGawiseret al.1 1200 ll 



HDFS 



FIG. 14. — Sky distribution of the 62968 objects in our photometric catalog 
(points) showing regions effectively masked out due to bright stars. 



which implies saturation in at least one filter. 

6. RESULTS 

6.1. R-band number counts 

The sky density of objects in our catalog is greatly re- 
duced near the brightest stars, leading us to create a mask 
to properly exclude these regions from analyses of sky den- 
sity and clustering. Our methodology is described in detail 
in L. Infante et al. (2005, in preparation). A careful analy- 
sis of the implied stellar and galaxy luminosity functions as 
a function of stellarity cut shows that stellar contamination is 

8 The flags column in the catalog is the result of a maximum performed 
over the flag values output by SExtractor for each measurement image run in 
dual-image mode. 




R 

FIG. 15. — i?-band galaxy counts per unit magnitude per square degree with 
Poisson error bars and solid fit of \og [0 (N) = — 3.52 + 0.34R based on counts 
at total AB magnitudes 20 < R < 24. Incompleteness is minimal below total 
magnitude R = 24.5 and reaches 50% around R = 25.5. 
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FIG. 16. — U — V versus V — R colors in our filter system (AB magni- 
tudes) of a Lyman break galaxy template spectrum (Adelberger, priv. comm.) 
from re dshift to 3.4 (solid curve), the Pickles catalog of stars (plusses, 
IPicklesll 998), and the definite point sources in our catal og (stellarity > 0.95 , 
points). The C and MC color selection regions used by Steidel et al. 12003) 
transformed into our UVR filter system are shown with dashed lines, and our 
adopted color selection region is shown with solid lines. 



iProchaska et al.ll2002l) . Figure [TBI shows star colors expected 
and measured in our survey versus our adopted UVR color 
selection criteria of 

(U-V)>1.2 

-1.0<(V-R) < 0.6 
(U-V)>3.8(V-R) + 1.2 

19<R < 25.5 , (13) 

where all magnitudes are AB. The color criteria shown as 
dashed lines in the figure result fr om the transf ormation of 
the U„ G 1Z criteria described by ISteidel et alJ i 19991) into 
UVR. 9 Although the upper branch of the stellar distribution 
at red V-R colors represents giants which are rare at our sur- 
vey depth and high galactic latitude, we shifted our criteria 
to avoid this region as it also contains dim dwarf stars with 
correspondingly large errors in color, which are the primary 
expected source of interlopers. The region of the Steidel et 
al. criteria avoided by this shift had the highest fraction of 
interlopers and also shows a high interl oper fraction in the 
more liberal t/V/?-selected subset of the iCooke et al.l (120051) 
LBG sample. The LBG template curve falls within our selec- 
tion region for 3.0 < z < 3.4. At higher redshifts, the LBGs 
rapidly become too red in V - R to distinguish from lower- 
redshift objects and are better selected as dropouts in B-R. 
We allow the color selection region to extend far to the blue 
in V - R to account for the effect of Lyman a emission lines, 
which add flux in V. 

For spectroscopic selection and clustering analysis, objects 
with U fluxes, fy, less than their ler flux uncertainty, ay, were 

9 For power-law spectra, the effective wavelengths of these filters imply 
the translations (U-V) AB = 1.25(I/„ -G) and (V-R) AB = 0.50(G-7£). 



assigned an upper limit of fy = ay to make it less likely that 
negative sky fluctuations were responsible for turning a dim 
object in U into a formal dropout in U — V. This correction 
avoids lower-redshift interlopers at the cost of some incom- 
pleteness. The vast majority of our LBG candidates are for- 
mal dropouts in the U filter, with fy < ay before being set to 
these limit values. The limit values were used to generate the 
U — V colors plotted in Fig. [n] which shows our full catalog 
including z — 3 Lyman break galaxy candidates selected by 
their UVR colors. 

Our selection criteria yield 1607 candidates i n 1137 square 
arcminutes, or 1.4 arcmin" ~ 2 . ISteidel et.aU {l999) found 
1.2 arcmin" 2 at 71 < 25.5 using stricter criteria than ours, 
ISteidel et alJ (|2003) found 1.7 arcmin 2 using the extended 
U„ G 7Z cri t eria sh own translated into UVR in Figure [TBI 
Cana k etafl J2j[)04) found 1.5 arcmin" 2 using UBR, and 
Cook e etaU (l2005|) fou nd 1.5 arcmin" 2 using u'VRI crite- 
ria. [Steidel et al. (2003) found a redshift distribution < z >= 
2.96 ± 0.29. Our UVR criteria are expected to start selecting 
LBGs at redshifts higher by 0. 16 because the red cutoff of the 
U -band filter and its effective wavelength after accounting for 
the MOSAIC II CCD QE are both ~ 200 A redder than for U„. 
Our criteria are expected to stop selecting LBGs at redshifts 
higher by up to 0.4 because the Lyman alpha forest causes red 
V-R colors at a higher redshift than in G-1Z due to the longer 
effective wavelength of V. Hence we expect < z >= 3. 2 ±0.4, 
and this will be calibrated via spectroscopy. 

Figure^Dshows the spatial distribution of these LBG can- 
didates. We esti mate the angular c orrelation function of these 
objects using the lLandv & Szalavl (IT991 estimator 



w(9) = 



DD-2DR+RR 



RR 



(14) 



where DD represents the number of pairs at separation 9 in the 
object catalog, RR represents the appropriately-normalized 
number of pairs at separation 9 in a Poisson-distributed mock 
catalog, and DR represents the appropriately-normalized 
number of pairs at separation 9 with one in the object catalog 
and the other in the mock catalog. This estimator has variance 



var[w(9)] 



d + w(g)) 2 
RR 



(15) 



where the denominator gives the Poisson error contribution 
i.e. the number of pairs expected to be found in the bin cen- 
tered on 9 given the survey geometry and the total number of 
LBGs in the survey. 10 The observed values of w(9) must then 
be corrected for the integral constraint caused by estimating 
the tru e sky density o f LBGs from our survey to find a power- 
law fit ( Fouca ud et al.l2003l) . 



A w 9-P=w ohs (9) + 



1 

O 2 



(16) 



In practice, this correction is performed by using a modified 
estimator 

(1 +A)DD-2DR + RR 

W( )= RR ' ( } 

where the correction factor A is given by 

J2H9)RR 



A = 



(18) 



A is estimated iteratively and usually converges in a few iter- 
ations. 

10 RR is equivalent to 1/p from Landv & Szalav 1993 
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FIG. 17. — U — V versus V — R colors of our full catalog (points) in AB magnitudes from SExtractor AUTO (Kron elliptical aperture) photometry (left panel) 
and from APCORR (corrected aperture) photometry (right panel). The LBG color selection region is shown with solid lines. The solid curve shows the colors of 
a Lyman break galaxy template spectrum over the interval < z < 3.4, with the interval 3.0 < z < 3.4 falling in the selection region. The reduced scatter reflects 
improved performance of the APCORR fluxes, with a clearer definition of the main sequence and M dwarf locus. 



The angular correlation function found for Lyman break 
galaxy candidates in Figure ^] is w(8) = (2.3 ± 1.0)(?~ ' 8 , 
where 9 is measured in arcsec, intermediate between previous 
results for fixed = 0.8 fr om Giav alisco et all (119981) of A w = 
1.4+^ and lFoucaud et alJ J2003I) of A,„ = 5.2 ± 0.7. This im- 
plies a correlation length for our sample of ro = 5 ± l/z~.'Mpc 
and an LBG bias factor of b ~ 2 for the ACDM cosmology. 
These results will be updated after spectroscopic calibration 
of our UVR LBG selection and by including our full square 
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FIG. 19. — Angular auto-correlation function of UVR-selected LBG can- 
didates (open circles with error bars) with fit (solid line) labeled. 
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FIG. 18. — Sky distribution of the 1607 LBG candidates in our catalog of 
E-HDFS showing some obvious clustering on small scales. 



7. CONCLUSIONS 

MUSYC is unique for its combination of depth and to- 
tal area, for the coverage in X-ray, UV, mid- and far- 
infrared wavelengths and for providing the UBVRIz' plus 
near-infrared photometry needed to produce high-quality pho- 
tometric redshifts over a square-degree of sky. This multi- 
wavelength coverage will enable comparison of selection ef- 
fects which have previously complicated the study of galax- 
ies in the high-redshift universe. The large area ensures good 
statistics for studies of clustering, luminosity functions, and 
surface densities. 

We expect that the careful attention paid to optimized pho- 
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tometry methods and the estimation of accurate photometric 
uncertainties will allow our public data to make a significant 
impact for a wide range of scientific studies. In particular, 
this high-quality UBVRIz' data, combined with forthcoming 
near-infrared photometry, should yield photometric redshifts 
with good accuracy. It is interesting to note that some of the 
standard techniques perform well despite making a number of 
idealized assumptions. For instance, we found that the opti- 
mal point source apertures are remarkably close in diameter 
to the idealized case of 1.35 fwhm due to a partial cancel- 
lation of the non-Gaussian PSF preferring a larger aperture 
and the correlated noise preferring a smaller aperture. The 
more traditional approach of using 2" apertures produces bi- 
ases in point-source colors. Alleviating those biases by "PSF- 
matching" i.e. smoothing the images to match the fwhm of the 
image with the worst seeing must be done carefully given the 
non-Gaussian PSF shapes and will result in degraded signal- 
to-noise for point sources. Estimates of point source detection 
depth that assume uncorrected noise and extrapolate from the 
rms pixel noise will overestimate the depth by about 0.5 mag- 
nitudes for 2" or 3" apertures. 

The APCORR fluxes introduced in this paper perform ex- 
tremely well, as seen in Fig. The nearly unbiased per- 
formance of APCORR fluxes for both point sources and ex- 
tended objects is impressive, and the correlated errors be- 
tween filters for extended objects are acceptable if the flux 
is misestimated by a constant fraction in all filters. How- 
ever, one should use caution when calculating photometric 
redshifts on extended objects in photometric catalogs gener- 
ated using apertures which might miss a different fraction of 



the object flux in each filter. The best results should be ob- 
tained by using the APCORR fluxes in our catalog for unre- 
solved and barely resolved sources, with further investigation 
needed to determine whether APCORR outperforms AUTO 
for extended sources. 

Our UBVRIz' images and catalogs of the 33' x 34.5' Ex- 
tended Hubble Deep Field South are available to the public at 
|http : / /www . astro . yale . edu/MUSYC| . 

We thank Buell Januzzi and Frank Valdes for their de- 
tailed answers and trouble-shooting with questions on the 
MSCRED package and Benne W. Holwerda whose guide, 
"Source Extractor for Dummies," provided a valuable sup- 
plement to the official SExtractor User's Manual written by 
E. Bertin. We acknowledge helpful conversations with An- 
ton Koekemoer, David Spergel, and Nick Suntzeff. We thank 
the additional members of the MUSYC collaboration, includ- 
ing Danilo Marchesini, Ken Rines, Ned Taylor, and Bill van 
Altena for their valuable input. We are grateful for support 
from Fundacion Andes, the FONDAP Centra de Astroffsica, 
and the Yale Perspectives in Science program. This material 
is based upon work supported by the National Science Foun- 
dation under Grant. No. AST-0201667, an NSF Astronomy 
and Astrophysics Postdoctoral Fellowship (AAPF) awarded 
to E.G. S.T acknowledges the support of the Danish Natural 
Research Council. We thank the staff of Cerro Tololo Inter- 
American Observatory for their invaluable assistance with our 
observations. This research has made use of NASA's Astro- 
physics Data System. 



APPENDIX 

OPTIMAL WEIGHTS FOR IMAGE STACKING FOR POINT SOURCES AND EXTENDED SOURCES 

We wish to combine a series of processed individual images which have signal Sj and noise N,. The exact signal obviously 
depends on which star is measured, but if a common set of stars are considered the variations in signal will be due to differences 
in exposure time, cirrus extinction, and atmospheric extinction. The noise should be sky-dominated away from bright objects but 
measuring the rms in "blank" regions of sky properly accounts for the readnoise as well. For surface-brightness-optimization we 
want to average the total signal from a set of stars, as measured by the IRAF routine phot or the inverse of the mscscale header 
output by mscimatch. For point-source-optimization we want to determine the average signal that falls within our eventual 
optimal aperture, which will have a diameter of roughly 1.35 times the fwhm of the final image. The seeing in the final image can 
be estimated or measured from an initially unweighted stack, or an iterative procedure can be used. Alternatively, a correction to 
the total star counts S'"' can be calculated by integrating an assumed Gaussian PSF out to the diameter of the expected optimal 
aperture of radius R; if the seeing in each image is a, = fwhm,/2.35, we obtain 

S^Sf'il-^-R 2 /!^)] . (Al) 

In either type of optimization, the noise value used should be the background rms for the expected aperture size, which depends 
on the level of noise correlation between pixels. Assuming that the noise correlations have the same behavior in each individual 
image, it is sufficient to use the pixel-by-pixel rms i.e. a\ measured in each image since this will produce the right relative 
weights. 

One can derive the weights that optimize the S/N of the combined image by noting that the combined image will have S = 
J2i w i$i an d N 2 = J2i w fN? (the weighted sum is just an unweighted sum of new images having signal vt/;Sj and noise W;iV;). 
Maximizing the function S/N produces the result w, = kSj/Nf, where multiplying all the weights by a constant k preserves the 
final S/N. In astronomy, it is typical to first scale the images to have equal signal levels; this accounts for differences in exposure 
time and extinction to provide constant photometry and makes it easier to delete outlying pixels due to cosmic rays, satellite 
trails, etc. The IRAF-mscred routine mscstack performs this scaling by multiplying by the mscscale header value and then uses 
a weightfile to weight the scaled images. So the weights used must be for the scaled images, which leads to the formula for 
surface-brightness-optimized weighting of 



1 (mscscalei x rms,) 2 ' 

which is equivalent to (Si/N/) 2 in the original image because the image being weighted is different than the one for which S, and 
N were measured. The output image is the same since instead of multiplying Si by Sj/Nf and adding we are equating the signal 
levels and then multiplying 1 by (Si/Nj) 2 , but performing the scaling first makes it easier to remove outlying pixels from this 
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sum as mentioned above. The weights being used correspond to S/N 2 for the scaled images; they reduce to the familiar case of 
inverse-variance weighting since the signal levels have already been equalized. For point-source optimized weighting we use 

wf s =( faC !° n V , (A3) 



v mscscale, x rmsi 
where factor-, is given by 

fact on = 1 - exp (-1.3 fwf ™«f ) , (A4) 
V fwhmf J 

where we have re-expressed the formula for 5, from Eq. lAll in terms of the measured seeing in an unweighted (or surface- 
brightness optimized) stack, fwhm stac k, and that in each individual image /w/iot,-. Again, these weights are only valid when used 
on the post-scaled images, and they differ by a factor of Sj from the weights appropriate for the original images. 

Due to the large number of individual exposures taken in e.g. U -band and multiple observing runs, it is sometimes necessary 
to perform stacking as an iterative procedure. To do this, the signal, noise, and seeing in the intermediate stacked images should 
be measured empirically as before. Then mscstack or the equivalent can be used to stack these intermediate images just as if they 
were individual images, scaling by the new mscscale values and calculating a new weights file using the above equations. As 
long as the final seeing is estimated well, no loss of signal-to-noise should occur from the iteration. 

A corollary question is when to discard an image with poor S/N rather than to include it in the stack. The formal answer 
appears to be never as w, > for all i. The magic of the optimal weighting formula can be seen: 



N 



IJ2S}/N? . (A5) 



Hence the optimal weights cause S/N to add in quadrature so it will never formally decrease no matter how poor an input image 
is. Note that having a magnitude of extinction or twice the seeing still leaves an image with significant S/N; a S/N threshold 
of 40% would require discarding images with about 2 times worse seeing or 1 magnitude of extinction (if the sky background 
is unaffected). The S/N of the combination of the image with a given S/N and an image with 40% of that S/N is 8% higher 
than that of the better image alone so this is somewhat useful. A conservative approach would be to cut entirely any images with 
S/N < 10% that of the median image as a way of reducing systematic effects not accounted for by these idealized formulae. 

OPTIMAL WEIGHTS FOR POINT SOURCE DETECTION 

Erwinl(IT985l) showed using simulations that the best performance for point source detection was achieved by filtering the image 
with the PSF itself. This result is derived in the SExtractor manual using Fourier transforms. In practice due to the constraints 
of computation time one chooses to cut off the PSF at a finite radius, as is standard in the convolution kernels offered as part 
of the SExtractor package. The general case for an optimized weighted sum of pixel values S = ~}2jWjSj, where each pixel has 

signal Sj and noise Nj, can be derived in analogy to the result for adding images found above, yielding wj oc Sj/N 2 . The constant 
background noise for sky-dominated objects thus leads to a filter of the precise shape of the PSF. For objects significantly brighter 
than the sky, S = N 2 giving constant weight i.e. a tophat which gradually morphs into the roughly Gaussian shape of the PSF as 
you move far enough away from the object center for the sky to dominate the noise. SExtractor only allows for a fixed convolution 
kernel so we optimized the detection for sky-dominated objects by using the PSF (truncated to 7x7 pixels) as our filter. 

By the same argument, the optimum signal-to-noise measurement of the photometry of a point source will be obtained by 
a weighted sum of pixel fluxes weighted using the PSF shape itself centered at the barycenter of the object. This is typically 
referred to as "PSF photometry" but it is not a supported feature of SExtractor photometry, just of object detection. Note that this 
is no longer strictly optimal when the noise is correlated between pixels as we have found in our images, since the derivation of 
the optimal weights for each pixel assumed uncorrected noise. We decided against using PSF photometry for this reason, and 
also because many of our science objects are slightly resolved and this would overweight the fluxes in their core, increase the risk 
of color biases, and complicate the task of correcting the aperture fluxes of extended objects. 
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TABLE 1 
MUSYC Fields 



Field 


RA 


Dec 


Galactic 


Ecliptic 


E(B-V) 


100/im Em. 


N(HI) 








Coords, (deg) 


Coords, (deg) 




(MJy/sr) 


(cm" 2 ) 


E-HDFS 


22:32:35.6 


-60:47:12 


(328,-49) 


(311,-47) 


0.03 


1.37 


1.6E+20 


E-CDFS 


03:32:29.0 


-27:48:47 


(224,-54) 


(41,-45) 


0.01 


0.40 


9.0E+19 


SDSS 1030+05 


10:30:27.1 


05:24:55 


(239,50) 


(157,-4) 


0.02 


1.01 


2.3E+20 


CW1255+01 


12:55:40 


01:07:00 


(306,64) 


(192,7) 


0.02 


0.81 


1.6E+20 



TABLE 2 

Optical Observations of E-HDFS with 
CTI04M+MOSAIC-II 



Dates 


Filter 


Number of 


Exposure 


Seeing 






Exposures 


Time (s) 


(") 


2002 Oct 6,8 


U 


47 


28200 


1.40 


2002 Oct 6,8,10 


B 


13 


7500 


1.35 


2002 Oct 6,10,12 


V 


30 


10440 


0.90 


2002 Oct 6,10,12 


R 


21 


6300 


0.85 


2002 Oct 6,12 


I 


26 


6300 


0.85 


2002 Oct 6,12 


z' 


11 


2700 


0.90 


2003 May 26,27 


U 


25 


15000 


1.20 


2003 May 27 


B 


10 


6000 


1.10 


2003 May 27,28 


R 


20 


6000 


1.10 


2003 May 26 


z' 


15 


3600 


1.30 
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TABLE 3 

Results of photometric calibration for 2002 October 6 in AB magnitudes 



Filter 


Airmass 


Seeing 
(") 


Zeropoint 


Airmass coeff. 


Color term 




U 


1.34 


2.1 


24.15±0.012 


0.55 


(-0.08±0.01)x(J/- 


-B)ab 


B 


1.38 


1.8 


25.24±0.009 


0.24 


(-0.03±0.004)x({/- 


V)ab 


V 


1.40 


1.6 


25.64±0.004 


0.15 


(0.03±0.003)x(fi- 


-R)ab 


R 


1.42 


1.5 


25.95±0.005 


0.09 


(-0.009±0.005)x(V 


-/)ab 


I 


1.44 


1.5 


25.40±0.006 


0.05 


(0.006±0.017)x(i? 


-/)ab 


z' 


1.46 


1.5 


24.61±0.014 


0.03 


(-0.06±0.054)x(fi 


-/)ab 



TABLE 4 

Total Exposure Time and Photometric Zeropoints for 
E-HDFS Images 3 



Filter Exposure Mag. Zeropoint Flux Zeropoint 

Time (s) (AB mag @ 1 count s~') (p]y (count/s)~' ) 



BVR 35340 24.29 0.698 

U 42600 23.35 1.660 

B 13200 24.85 0.417 

V 10440 25.40 0.251 

R 11700 25.75 0.182 

I 6000 25.27 0.283 

z' 6000 24.39 0.637 



"before correction for E(B-V)=0.03 



TABLE 5 
Fit parameters for 
background fluctuations 
in counts i" 1 as a function 
of aperture size using eq. 

El 



Filter 


ft 


a 


P 


BVR 


0.0029 


0.68 


1.47 


U 


0.0012 


0.82 


1.30 


B 


0.0043 


0.88 


1.33 


V 


0.014 


0.69 


1.34 


R 


0.019 


0.69 


1.40 


I 


0.032 


0.72 


1.42 


z' 


0.030 


0.87 


1.41 



TABLE 6 

Seeing, Apertures, and Source Detection Limits for E-HDFS Images 



Filter 


IRAF fwhm 
Mode (") 


SE fwhm 
Median (") 


SE r l/2 
Median (") 


Optimal Aperture 

(") 


Flux Enclosed 
(Fractional) 


Source Detection 
Limit (5 <r, AB) 


BVR 


0.95 


0.99 


0.59 


1.2 


0.48 


26.3 


U 


1.30 


1.48 


0.78 


1.6 


0.50 


26.0 


B 


1.25 


1.29 


0.71 


1.4 


0.47 


26.1 


V 


0.90 


0.96 


0.56 


1.2 


0.52 


26.0 


R 


0.90 


0.97 


0.58 


1.2 


0.49 


25.8 


I 


0.80 


0.99 


0.56 


1.2 


0.49 


24.7 


z' 


0.90 


1.06 


0.62 


1.2 


0.43 


23.6 



MUSYC: Images and Catalogs of E-HDFS 



TABLE 7 

Optical source catalog: astrometric and photometric parameters determined from BVR detection image and UBVRh! 

PHOTOMETRY. 



Number 


Name 


KA 


Dec 




y 


Mellanty 


r l/2 


A_l_IVLALrE, 




MUSYC- 


(hrs) 


(deg) 








( ) 


( ) 


B_IMAGE 


y_IMAGE 


FLAGS 


/_auto_U 


<r_auto_U 


/_auto_B 


cr_auto_B 


/_auto_V 


tr_auto_V 


/ti\ 


ccw from X-axis 




(tdy) 


OJy) 


0*Jy) 


0*Jy) 


(MJy) 


(/iJy) 


/_auto_R 


<r_auto_R 


/_auto_I 


cr_auto_I 


/_auto_z' 


<7_auto_z' 


/_apcorr_U 


cr_apcorr_U 


/"_apcorr_B 


0*Jy) 


(/*Jy) 


(MJy) 


(AtJy) 


0*Jy) 


OJy) 


(MJy) 


(MJy) 


(MJy) 


cr_apcorr_B 


/_apcorr_V 


a_apcorr_V 


/_apcorr_R 


cr_apcorr_R 


/_apcorr_I 


<r_apcorr_I 


/_apcorr_z' 


<T_apcorr_z' 


OJy) 


OJy) 


(pJy) 




0*Jy) 


0*Jy) 


(MJy) 


(pJy) 


(MJy) 


o 


T773073 68-60481 7 4 


97 506580 


-60 8048^0 


7397 04 


361 3 S3 


4Q0000 


377004 


0771 630 


o fi77 i 630 


4S 0000 


16 


00SQQQQ8 


01 73^Q 


0768763 

U.UZ.UOZ.UJ 


01 33876 


0771 88Q 


01 8Q7^3 


f) 053571 5 


01 92685 


0584936 


0536170 


1 05203 


0.142492 


01 531 94 


0299205 


0582284 


o 0300677 

U .U JUV7U / Z. 


1 99481 


0329477 


0QQ4549 


0346033 


0063358 

U.U7UJJJ O 


0964758 


129852 


999986 
V7 . z y zz. o u 


\ 


T993094 98-604008 Q 


99 506747 

Z,Z..JUU / 1 / 


-60 6601 50 


7395 69 


5443 56 


380000 


381 543 

Lf . JO IJtj 


0771 630 


0771630 


45 0000 


16 


0545032 


0134836 


0528900 


0133148 

U.UJJJ17U 


0600472 


01 81627 

vj.vj ioi uz / 


0531 777 

UivJJ L 1 I 1 


01 89471 


0346594 


053731 7 


-0 0605 1 77 


n 0999073 


173143 


0319547 


0.111140 


0292123 


0983857 


0301251 


0985063 


0331430 


0733320 

\J,\J 1 ~J~J^f^\J 


0938256 


-0 113113 


287106 


2 


J223023.85-604558.3 


22.506627 


-60.766210 


7397.14 


4134.88 


0.490000 


0.555360 


0.216003 


0.133500 


-0.200000 


24 


0.0909909 


0.0243018 


0.0549541 


0.0255182 


0.0440691 


0.0353969 


0.130563 


0.0386372 


0.159939 


0.108677 


0.0605380 


0.286208 


0.100151 


0.0285476 


0.0699377 


0.0270144 


0.0861040 


0.0280836 


0.125653 


0.0311628 


0.113170 


0.0863858 


0.0282574 


0.261025 


3 


J223024.98-603027.9 


22.506940 


-60.507770 


7394.68 


7619.57 


0.940000 


0.328410 


0.199716 


0.120951 


18.8000 


24 


-0.0217973 


0.0159112 


0.000543353 


0.0123818 


0.115021 


0.0169463 


300.840 


0.0511371 


577.641 


0.128364 


-0.00643368 


0.124763 


-0.0458961 


0.0422158 


-0.000965204 


0.0289711 


0.222445 


0.0300914 


601.662 


0.125434 


1191.94 


0.301193 


0.00430643 


0.280111 


4 


J223023.35-605242.2 


22.506487 


-60.878400 


7397.96 


2622.20 


0.610000 


0.613032 


0.369528 


0.202653 


-5.40000 


24 


-0.00917560 


0.0288420 


0.0469781 


0.0383380 


0.134234 


0.0542395 


0.155328 


0.0622344 


-0.104166 


0.130648 


0.110060 


0.458585 


-0.0224367 


0.0274903 


0.100572 


0.0275703 


0.155969 


0.0290837 


0.143652 


0.0332477 


0.0475046 


0.0869854 


0.677608 


0.276257 



NOTE. — This table is published in ils entirety in die electronic edition of the Aslrophysical Journal. A portion is shown here for guidance regarding its form and eontenl. Due lo geometry, ihese iirsl few objects are near the image border 
and have poor photometry, leading to their large FLAG values. 



